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ABSTRACT 

We have recently measured the angle between the spin and orbital angular momenta of PSR 
B1534+12 to be either 25 ± 4° or 155 ± 4°. This misalignment was almost certainly caused by 
an asymmetry in the supernova explosion that formed its companion neutron star. Here we combine 
the misalignment measurement with measurements of the pulsar and companion masses, the orbital 
elements, proper motion, and interstellar scintillation. We show that the orbit of the binary in the 
Galaxy is inconsistent with a velocity kick large enough to produce a nearly antialigned spin axis, so 
the true misalignment must be ~ 25°. Similar arguments lead to bounds on the mass of the com- 
panion star immediately before its supernova: 3 ± IMq. The result is a coherent scenario for the 
formation of the observed binary. After the first supernova explosion, the neutron star that would 
eventually become the observed pulsar was in a Be/X-ray type binary system with a companion of 
at least 10-12 Mq. During hydrogen (or possibly helium) shell burning, mass transfer occurred in a 
common envelope phase, leaving the neutron star in a roughly half-day orbit with a helium star with 
mass above ~ 3.3M Q . A second phase of mass transfer was then initiated by Roche lobe overflow 
during shell helium burning, further reducing both the helium star mass and orbital period before the 
second supernova. Scenarios that avoid Roche lobe overflow by the helium star require larger helium 
star masses and predict space velocities inconsistent with our measurements. The companion neutron 
star experienced a velocity kick of 230 ± 60 km/s at birth, leading to a systemic kick to the binary of 
180 ± 60 km/s. The direction of the kick was roughly opposed to the instantaneous orbital velocity 
of the companion, but the kick angle is largely unconstrained. 
Subject headings: pulsars: individual (PSR B1534+12) — supernovae 



1. INTRODUCTION 

The handful of known short-period double neutron star 
binaries are of interest to physicists and astrophysicists 
working in a wide variety of fields. Most famously, stud- 
ies of two of these systems have been used as labora- 
tories to test aspects of fundamental gravitation theory 
that are inac cessible in the solar system, including the 
prediction bv lEinst emjjLjnl ^ of the production of gravi- 
tation al radiation (Tavlor & Wc isberglll98 9: Stair s et al.l 
|2002|) . Such binaries, in their late stages of gravita- 
tional inspiral, are among the most likely sources to 
be detected by LIGO and othe r gravitational wave ob- 
servatories ({Abbott et al]l2004l and references therein). 
They yield high-precision mass measurements that have 
been used to study the neutro n star equation of state 
l)Thorsett fc Chakrabartvl lT999). Finally, they are im- 
portant long-lived relics of binary evolution that pre- 
serve in their orbital elements a memory of short-lived 
mass transfer phases that are rarely if ever observed in 
progress (e ; g.. iBhattacharva fc van den Heuvell 19911). 

A complicating factor in understanding the birthrate 
of these binaries has been constraining asymmetries in 
supernova explosions. Isolated radio pulsars are ob- 
served to have high space velocities, interpreted as ev- 
idence that newborn neutron stars receive a momentum 
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kick at birth (jDewev fe Cordesl Il987t iLvne fe Lorimerl 
Il994|) . There is some evidenc e that kicks given 
to stars in binaries are smaller !G^irrjJ^_e^alJ |129_4|; 
iCordes fc ChernofB 119 971 iPortegies Zwart fc Yung clson 



1998; Hughes fc Baileslll999|) . but it is unclear whether 
this is a selection bias related to the probability of binary 
survival, or whether these stars truly have a different kick 
distribution, perhaps because of a different angular mo- 
mentum at collapse. Since kicks introduce new degrees of 
freedom in the change of the binary elements after the su- 
pernova event, understanding kicks is also important for 
estimating the pre-supernova orbital properties. Finally, 
asymmetries can produce misalignments between the 
spin and orbital angular momenta in the binary, which 
can lead t o changes in the gravitational-wave signature 
at m erger l)Apostolatos et allll994t iGrandclement et alJ 
2003). Misalignment also makes possible the study 
of general relativistic (geodetic) precession, providing 
a q ualitatively new test of strong field gravity theo- 
ries llDamour fc RuffinillT974) iBarker fc Q'Connell|T975| 
i Barker fc O'Connelll 119751: lEsoosito fc Harrisonl 119751: 
IDamour fc Tavlorlll992t IStairs et al.H2004n . 

Recently, we have used precession measurements (to- 
gether with the assumption that general relativity cor- 
rectly describes the kinematics of neutron star binaries) 
to meas ure the misalignme nt angle in the PSR B1534+12 
binary l)Stairs et al] 12001) . By combining this result 
with our previous high precision measurements of the or- 
bital elements, component masses, system distance, and 
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V 

To observer 

Fig. 1.— The geometry of the PSR B1534+12 orbit. The or- 
bital angular momentum vector is k, and the spin axis is S. The 
precession cone, with half-opening angle <5, is shown. The vector i 
is directed towards the ascending node. 



proper motion, and with a measurement of the orien- 
tation of the orbit from interstellar scintillation studies, 
we are able to tightly constrain the magnitude and di- 
rection of the asymmetric velocity kick imparted to the 
pulsar's companion when it was formed. Three unknown 
parameters — the line-of-sight velocity of the binary, the 
pre-supernova orbital size, and the pre-supernova mass 
of the companion — are all tightly constrained. These in 
turn strongly constrain the evolutionary history of the 
binary, at least during the period between the first and 
second supernovae, as well as the supernova asymmetry. 

Our work here draws on previous studies of several rel- 
ativistic binaries, particularly of this system and of PSR 
B1913+16. However, the quality and extent of the data 
for PSR B1534+12 now far exceeds what is available for 
the other binaries, allowing much more robust conclu- 
sions to be drawn. We describe these data in In 
we review the generic evolution of a compact binary sys- 
tem before, during, and after the second supernova event. 
In 21 we combine observations with theory to constrain 
the particular history of the PSR B1534+12 binary. 

2. OBSERVATIONS OF PSR B1534+12 

As we will see below, the current orbit and kinemat- 
ics of the PSR B1534+12 binary were largely determined 
by the properties of the binary just prior to the second 
supernova event (which created the neutron star com- 
panion to the observed pulsar) and by the mass loss and 
momentum imparted to the system during the second 
supernova. Neither the size nor the eccentricity of the 
orbit have changed significantly since that event, while 
the binary has moved ballistically through the Galac- 
tic gravitational potential and gyroscopically maintained 
the direction of its orbital angular momentum vector. A 
wealth of data fully constrain the current orbital prop- 



erties, including the orientation of the orbit on the sky. 
(See Fig. n for a reference to angles describing the orien- 
tation.) The distance and transverse components of the 
binary velocity are also well constrained. Only the ra- 
dial velocity and the age cannot be directly determined, 
though a rough upper limit can be placed on the time 
since the second supernova. 

2.1. Orbital elements and masses 

Radio pulsar timing has produced exquisitely accurate 
estimates of the orbital elements of the PSR B1534+12 
binary. The system is essentially a single-line spectro- 
scopic binary, with the Doppler shift of the pulsar pe- 
riod replacing the usual Doppler shift of a stellar spec- 
tral line. Although only line-of-sight velocity variations 
can be observed through timing, the measurement of 
general-relativistic corrections to the Keplerian orbital 
equations allow the orbital inclination and component 
masses to be determined. The timing experiments, anal- 
ys is, and resul t s hav e been described in extensive detail 
in lStairs et al.1 (|2002|) . For our current purposes, it is suf- 
ficient to know the orbital period P& = 0.421 days, eccen- 
tricity e = 0.274, relative semimajor axis a — 7.62 light 
seconds, orbital inclination sini = 0.975, pulsar mass 
mi = 1.333Mq, and companion mass mi — 1.345A/q. 
In all cases, the uncertainty is small relative to the last 
quoted digit. 

2.2. Misalignment of spin and orbital angular momenta 

Misalignment between the pulsar spin axis and the or- 
bital angular momentum axis causes the spin axis to pro- 
cess around the orbital axis at a rate of 0.51°/yr. This 
so-called geodetic precession is a simple relativistic con- 
sequence of the parallel transport of the spin vector in 
curved space; it is a consequence of the same angle deficit 
that produces precession of the periastron of the orbit. 

Precession of the pulsar spin axis causes the impact 
angle between the observer's line of sight and the pul- 
sar's magnetic pole to change, with an amplitude that 
depends on the misalignment angle S between the spin 
and orbital angular momenta. This shifting viewing an- 
gle causes changes in the apparent pulsar beam shape 
and polarization properties. 

Recently, we have used a long-term study of the po- 
larization and bea m shape to measure S — 25.0 ± 3.8° 
l|Stairs et al J 12004]) . There is a discrete ambiguity that 
corresponds to reversing the spin direction, so S = 
155.0 ± 3.8° is also allowed. 

2.3. Position and velocity 

In addition to the orbital elements, timing measure- 
ments provide very accurate estimates of the pulsar po- 
sition and proper motion, through the annual modula- 
tion of the observed pulse period caused by the Earth's 
motion. In the case of PSR B1534+12, the distance can 
also be estimated from a measurement of the apparent 
orbital period derivative (after correction for the contri- 
bution due to gravitational radiation emission) — which 
is caused by a combination of acceleration in the Galac- 
tic potential and apparent acceleration along the line of 
sight that arises from the transverse motion. 

The binary's current position is well out of the Galac- 
tic plane, at I = 19.8°, b = 48.3°. The distance is 
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d = 1.02 ± 0.05 kpc. The observed prop er motion is 
U£A = 1.3 mas/yr, fi dcc = -25.1 mas/yr l|Stairs et all 
2002). The total proper motion is /i = 25.2 mas/yr, di- 
rected towards a Galactic position angle 239.2°, which 
at the known distance implies a transverse velocity (un- 
corrected for Solar motion) of 122 km/s. In contrast to 
the usual case with a spectroscopic binary, there are no 
direct constraints available on the radial component of 
the space velocity as the rest-frame spin period of the 
pulsar is unknown. 

In projection, the pulsar appears to be moving towards 
the plane, with a ^-component velocity —62 km/s. How- 
ever, a radial velocity component of just 56 km/s is suf- 
ficient to reverse the apparent motion, so it is unknown 
whether the pulsar is presently moving towards or away 
from the plane. 

2.4. Scintillation 

Timing studies, which measure the component of ve- 
locity along the line sight, are insensitive to rotations 
of the binary orbit around the line of sight, and hence 
the angle f2 of the ascending node on the plane of the 
sky. 1 This angle is accessible through proper- motion- 
induced changes in the projected orbital size, with 
x/x = /isinjcoti, with i the orbital inclination and j 
the angle between the pr oper motion direction and Q 
ijArzoumanian et al.lll996j) . Because fi is small for PSR 
B1534+12, this has not yet been measurable. 

In principle, some constraint on the direction of the 
pulsar's spin axis could be obtained from polarization 
studies, if corrections are made for Faraday rotation and 
observations are properly calibrated against an absolute 
reference source. This has also not yet been done for 
PSRB1534+12. 

An alternative approach that has been successful for 
this source comes from study of diffractive scintillation. 
The propagation of the pulsar signal through the in- 
homogeneous interstellar medium produces a pattern 
( "screen" ) of intensity variations caused by construc- 
tive or destructive interference. The speed at which 
this screen sweeps past the observer is proportional 
to the vector sum of the Earth's motion, the binary's 
proper motion, and the pulsar's orbital motion, and is 
inversely proportional (at least statistically) to the decor- 
relation timescale for intensity variations ijLvnd 119841: 
IDewev et aflll98c| . 

By measuring the decorrelation timescale as a func- 
tion of orbital phase, iBoedanov et alJ l)2002l) were able 
to measure Q = 70 ± 20° or 290 ± 20°, reckoned north 
through east. The solutions correspond to cosi < and 
cosi > 0, respectively; 2 o nly the second is con sistent with 
our polarization studies ijStairs et alJl2004j) . Although 
the uncertainties remain large, we will see that the re- 
sulting constraint on the direction of the orbital angular 
momentum is useful. 

2.5. Age 

1 Note that we use the standard convention in binary pulsar 
studies and define the "ascending node" as the node at which the 
observed pulsa r is moving away from t he observer. 

2 Although IBogdanov et all 12002:1 do not explicitly identify 
which solution corresponds to which sign of cosi, our independent 
analysis has been confirmed by Bogdanov, private communication. 



A simple estimate for the age of the binary since the 
second supernova is given by the characteristic age of the 
observed pulsar: r c = P/2P = 0.25 Gyr. The true age 
r is unknown. If the pulsar was born at or below the 
spin-up line and has evolved by magnetic dipole brak- 
ing (with braking index n — 3), then r ~ 0.21 Gyr 
l)Arzoumanian et alJ Il999[) . The age could be signifi- 
cantly less than this if the pulsar spin-up ended before 
the spin-up line was reached, or if mass transfer ended 
(and spin-down began) well before the second supernova, 
and could be somewhat larger if n < 3. As will be dis- 
cussed in H4.2I ages as small as 10 Myr allowed if the 
binary is just leaving the plane for the first time. If the 
system is truly that young, the companion star may still 
be an active radio pulsar, which cannot be ruled out if the 
Earth lies outside the path of its lighthouse beam. We 
believe it is unlikely that the age of the binary is much 
larger than the characteristic age. The direct constraints 
on effective braking indices for pulsars with intermediate 
mag netic fields, lik e PSR B1534+12, is limited. How- 
ever, iBrisken et alJ l)2003|) have found that characteris- 
tic ages and kinematic ages agree well in a set of high 
field pulsars, and in a sa mple of millisecond p ulsars (with 
smaller magnetic fields') . ICamilo et alJ (JT994) found that 
half had characteristic ages larger than a Hubble time, 
implying that characteristic ages were typically an over- 
estimate of true age. Without further information, we 
take 2.1 x 10 8 yrs as a rough upper limit to the age of 
the system, but none of our conclusions depend strongly 
on that limit. 

3. THE EVOLUTION OF THE PSR B1534+12 BINARY 

In general terms, the formation of close double neu- 
tron starjDinarijj^yfice^^ understood 
Csee iBhattacharva fc van den Heuvellll99lL for example, 
for a review). It begins with a first supernova in a 
binary star system, in which the neutron star that is 
now observed as a pulsar was born. After a possible 
Be/X-ray binary phase, a period of mass transfer occurs: 
"case B" if it occurs by Roche lobe overflow during hy- 
drogen shell burning, as the companion first climbs the 
giant branch, or "case C" if it occurs during helium shell 
burning. ("Case A" transfer, while the companion is 
still on the hydrogen main sequence, is thought to lead 
to destruction of the binary through merger.) During 
the ensuing common envelope evolution, the system spi- 
rals together, leaving a binary containing a neutron star 
and a bare helium star. After case B transfer, the he- 
lium star is unevolved. After case C transfer, the helium 
star has already developed a carbon-oxygen core. A sec- 
ond stage of mass transfer may occur, which is again 
called case A, B, or C if it occurs during helium core 
burning, helium shell burning, or carbon burning, respec- 
tively. The overall mass-transfer history is, for example, 
case BB if a first phase occurred during hydrogen shell 
burning and a second during helium shell burning. In 
any case, during mass transfer the orbit is circularized 
and the pulsar is spun-up, or "recycled," and the pulsar 
spin axis aligned with the orbital angular momentum. 
Finally, a second supernova explosion produces the com- 
panion neutron star, which may live for ~ 10 7 yrs as a 

pulsar. (The faster, recycled pulsar can live ~ 10 9 yrs, 
so most observed double-neutron star binaries are simi- 
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lar to B1534+12, with only the recycled pulsar visible, 
and only in rare cases like the recently discovered PSR 
J 0737-3 39 llLvne et al.l2004|) can both pulsars be seen.) 

iBrownl l)1995|) proposed an alternative scenario in 
which the initial inspiral of two stars of similar mass 
produces a double helium star binary. Though differ- 
ent in the early phases of evolution, this scenario still 
requires mass transfer after the first supernova in order 
to recycle the observed pulsar, so the system just prior 
to the second supernova is again expected to be a close, 
circular neutron-star-helium-star binary with the pulsar 
spin and orbital angular momenta aligned. We will not 
further discuss this possible evolutionary track, except 
to note that because the late evolution is very similar to 
the standard scenario most of our conclusions would be 
unchanged. 

At this point, it is important to emphasize that the 
immediate progenitors of systems like PSR B1534+12 
that is, close binaries that contain a neutron star and 
bare helium star — are not observed, with the possible 
exception of Cyg X-3, which appears to be a Wolf-Rayet 
star in orbit around either a neutron star or a black hole 
l|van Kerkwiik et alJll992T: iFender et~all fl9991. Our pri- 
mary observational constraints on the late stages of the 
standard evolutionary scenario are the properties of the 
resulting double neutron star binaries, including their 
masses, orbital elements, and space motions. In the 
case of PSR B1 534+12, a number of studies have al- 
ready been done ^ fmn^n^L^iL^ill22i\EasLlLKa^p ffeT4 
19971 lArzoumanian et a,l.lll999t iFrancischelli et a,l.ll2002t 
Dewi fc Pols! 120031 iWillems et al.ll2004|) . To these, we 
now add impro ved estimates of th e pulsar distance and 
proper motion ijStairs et al.ll2002[). measure ment of the 
spin-orbit misalignment l|Stairs et all 120041). and mea- 
surem ent of the line of nodes angle ( Bogd anov et al.l 
2002), which together strongly constrain the magnitude 
of any asymmetric kick applied to the companion star 
during the supernova in which it was born. 

3.1. The pre- supernova binary 

A number of constraints on the properties of the pre- 
supernova binary come from the size of the helium star 
companion, the core mass required to form a neutron 
star, and the possible stellar evolution and mass trans- 
fer histories that could produce the close binary. A de- 
tailed study of the evolution of close helium sta r bina- 
ries has recently been c arried out by iDewi et aD (|2002f) 
and lDewi fc Pols! ((2003), and we draw heavily from their 
work. 

The small size of the pre-supernova binary — far too 
small to contain two main sequence stars — immediately 
implies that the orbit went through at least one mass 
transfer and inspiral phase. In the standard model, the 
system after the first supernova is a Be/X-ray type bi- 
nary that spirals together in a common envelope phase 
during either case B or case C mass transfer. (As noted 
above, an alternate model posits inspiral during a double 
helium star phase before the first supernova.) 

The subsequent evolution depends on the helium star 
mass and whet her a second mass transfer phase occurs. 
As discussed in IDewi et alJ l)2002j) . case BA mass trans- 
fer, during helium core burning, never results in a double 
neutron star binary. Low mass stars lose enough mass 
that they leave CO white dwarf remnants, rather than 



collapsing to neutron stars. High mass stars experience 
dynamically unstable Roche lobe overflow, and probably 
coalesce with the neutron star to form an isolated black 
hole. 

Case BB mass transfer, during helium shell burning, 
results in the removal of the helium star's envelope and 
the formation of a CO or ONe white dwarf if the zero age 
helium star mass is less than 2.5M Q . Stars with mass 
less than 3.3M (or 3.8M in orbits below about six 
hour period) spiral in during a common envelope phase, 
and probably produce very compact (15 minute) double 
neutron star binaries, which have ver y short merger time s 
and are therefore difficult to detect ijDewi fc Polsl f2003). 
More massive stars, up to about 6.5M Q , avoid a common 
envelope during Roche lobe overflow, and hence avoid 
catastrophic inspiral. After mass transfer, these systems 
have orbital periods of hours to days, and pre-supernova 
core masses between about 2.2 and 4M Q . As we will see, 
these are progenitors for systems like PSR B1534+12. 

Larger helium stars don't swell significantly as they 
evolve, so avoid Roche lobe overflow. Binaries contain- 
ing such massive He stars never experience Roche lobe 
overflow, and are also possible progenitors of double neu- 
tron star binaries. In this case, the spin-up of the ob- 
served pulsar occurs by accretion from the helium star 
wind. In this scenario, a lower limit on the zero age he- 
lium star mass comes from the requirement that it never 
expand into contact with its Roche lobe, and an upper 
limit comes from the requirement that its zero age ra- 
dius fit within t he Roche lobe. In the specific case of 
PSR B1534+12. [Dewi fc Polsl (|2003|) have found that al- 
though a fairly wide range of zero age helium masses are 
possible (up to about 13M Q ), after wind losses there is 
only a very limited range of pre-supernova mass possible: 
5 ± O.5Af . As we will see, supernovae in these systems 
can reproduce the observed orbit of PSR B1534+12, but 
only with a kick too large to be consistent with the ob- 
served space velocity. There are no possible scenarios for 
the formation of PSR B1534+12 that begin with case C 
transfer and then avoid Roche lobe overflow. 

3.2. The second supernova event 

We denote the orbital eccentricity, relative semimajor 
axis, and total mass before the second supernova by ej, 
ai, and Mi — mi + m2j, where m, is the mass of the 
observed pulsar and m^i is the presupernova mass of the 
companion. Before the second supernova explosion, the 
mass transfer that spun-up PSR B1534+12 also almost 
certainly circularized the binary orbit, so we take ej = 0. 
The orbital velocity just prior to the supernova is "V%\ 
from Kepler's laws the (constant) orbital speed is given 
by = GMijai. During the explosion, mass M; — 
Mf = mn — mi f is lost from the system, and a velocity 

kick Vfe (the "natal kick") is applied to the companion 
star. After the explosion, the new elements are e/, af, 
and Mf = mi + m2f, and the change in the direction 
of the orbital angular momentum is S. E xpress i ons fo r 
the postexplosion elements can befoundin[HnilCl983). 
For convenience, we follow the very clear exposition of 
iWex et all l)2000() and define normalized quantities v = 

Vk/Vi, a ee df/di, = Mf/Mi, and r\ = ^ a(3{\ - ej). 

Then the shifts in the elements can be related to the 
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components of the kick velocity vector through: 



and 



■■±r, 



1 



"(! + e /)/ \ a ( l ~ e /) 



1/2 



v y = t] cos 5—1 
v z = ±r/ sin S. 



(1) 

(2) 
(3) 



Here the x axis is the line connecting the pulsar and 
companion at the time of the explosion, y is the direc- 
tion of instantaneous velocity of the companion, and z is 
perpendicular to the orbit. From equation^it is straight- 
forward to see that the initial semimajor axis is bounded 

by 

1 1 

< a < . (4) 



1 + ej 



1 



"f 



For any a in this range, j3 < 1, and S, we can use equa- 
tions. 1° calculate the kick velocity vector needed to 
produce the observed system. 

Given the mass loss and kick applied to the compan- 
ion star, the resulting kick to the binary center of mass 
u = u/Vi can also be calculated. We refer to this as the 
"systemic kick." (Note that the binary receives a sys- 
temic kick even i n a symmetric ex plosi on, when V k = 0.) 
Again, we follow iKalogeral l|1996ft and lWex etafl l|2000li 
and write 



u 2 = Ki + k 2 (2 — a 1 ) — K3 



^Ja ^1 - e'jj cos S, (5) 



where 



_ ,u 2f 



= — ^, no — : — , and «q = 2 a /ki ko- (6) 

M?' M l M f ' d V V ' 

We can decompose this into the component perpendic- 
ular to the post-supernova orbit u±_ and the component 
in the plane of the orbit fin, with 



Ki sin 6 



(7) 



and 



u\ = Ki COS 2 



5 + k 2 (2 — a x ) — K3^Ja ^1 — e 2 ^j cos 5, 

(8) 



3.3. Evolution after the second supernova 

The current binary system consists of two neutron 
stars. The orbit evolves slowly as it loses energy and 
angular momentum to gravitational radiation. Neither 
mass transfer nor tidal effects will have any significant 
contribution to the orbital evolution until the last stages 
of inspiral, about 3 Gyr from now. 

Although high precision radio pulsar timing observa- 
tions have measured the effects of gravitational radiation 
damping on the orbital elements, they are very small. In- 
deed, the total change in the elements since the second 
supernova explosion formed the double-neutron star bi- 
nary c an for most purposes be ignored. Following iPetersI 
( 1964), we can write the fractional rates of change of the 
semi-major axis and eccentricity as: 



Ida _ MG A m 1 m 2 M 

a~dt ~~ 5C 5 a 4(!_ e 2)7/2 



73 



37 



— -I e 

24 96 



(9) 



1 de 
e dt 



304G 3 mi m 2 M 
' 15c 5 a 4(i_ e 2) 5 / 2 



121 

304' 



(10) 



For PSR B1534+12, the timescales for evolution of a 
and e are 9.0 Gyr and 2.0 Gyr, respectively: both much 
longer than the pulsar characteristic age. We can there- 
fore ignore post-supernova evolution of the orbital ele- 
ments. As we will see, we are actually interested only in 
a and in lie. Although we can't properly correct these 
quantities to their original values without knowing the 
true binary age, the errors that enter when we neglect 
radiation damping are ~ 3%, which will be dwarfed by 
uncertainties in stellar modeling. 

Because tidal effects arc extremely small, no 
gravity-driven evoluti on of 5 is expected. Recently, 
IDemorest et al.l l|2004fl have suggested that wind from 
the fast pulsar in the PSR J0737— 3039 binary, which 
penetrates well into the light-cylinder of the slower pul- 
sar, may have aligned the spin axis of the slow pulsar 
with the orbital angular momentum. Such alignment is 
unlikely for a fast pulsar like B1534+12; for any likely 
birth luminosity, the standoff radius would have been 
well outside the light cylinder of B1534+12. The orbital 
motion dominates the total angular momentum of the 
system, so the orientation of the orbit remains nearly 
fixed in space, while the spin axis of the pulsar precesses 
around it tracing a cone with half opening angle 6. 

In summary, because of its small age and slow evo- 
lution, we can without significant error take the cur- 
rently observed orbital parameters, orbital orientation, 
and spin-orbit misalignment to be an accurate record of 
the state of the system immediately after the second su- 
pernova explosion. The space velocity imparted to the 
system during the second supernova is not, however, pre- 
served as the binary moves through the Galactic poten- 
tial. We discuss this complication in H4.2I 

4. DISCUSSION 

Our goal is to combine the wealth of observational data 
on the current state of the binary system and its po- 
sition and motion through the Galaxy with constraints 
that come from binary evolution modeling and the effects 
of the mass loss during the supernova event to develop 
a consistent model of the late evolution of the system. 
Several free parameters remain that must be constrained. 
Among the quantities not directly accessible from obser- 
vations are the radial velocity (and hence the full space 
velocity), the age and birthplace, the pre-supernova mass 
and radius of the binary orbit, and the asymmetry of the 
supernova explosion. The constraints on these parame- 
ters are interrelated and complex. 

To clarify the situation, we begin by allowing the as- 
sumed radial velocity of the system to vary. We will find 
that assuming the system was born near the plane of the 
Galaxy (where most massive stars are found) in the last 
210 Myr is sufficient to severely limit both the current 
space velocity of the pulsar and the kick velocity given to 
the binary system at birth. This systemic kick, the post- 
supernova eccentricity and size, and the misalignment all 
depend on the pre-supernova orbital size, the mass loss 
during the supernova, and the magnitude and direction 
of any asymmetric kick given to the companion. 
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We will show that the range of kick magnitudes and 
progenitor masses that satisfy all the constraints is very 
limited. 

4.1. The binary elements 

The relatively low eccentricity of the observed bi- 
nary together with equation 0] constrains the size of 
the pre-supernova binary to within about 30%. Evo- 
lutionary arguments constrain the pre-supernova mass 
of the companion star to between about 2 and 5.5 solar 
masses. Precession studies constrain the post-explosion 
misalignment with the discrete a ambiguity between 
nearly aligned pulsar spin and orbital angular momenta 
(5 ~ 25°) and nearly counteraligned angular momenta 
((5-155°). 

Given the observed orbital parameters for 
PSR B1534+12, we can rewrite equations ^ HO as 

V x = ±380 [{a - 0.785)(1.377- a)] 1/2 (11) 

V y = 395 (±0.872a - y^j (12) 

V a =±160a (13) 

In eauation ll2l the + sign refers to the nearly aligned 
case, the — sign to the nearly counteraligned case. 

From the observed parameters, we know that 0.785 < 
a < 1.38. With no assumptions about the mass loss 
in the second supernova explosion, this gives us di- 
rect constraints on two components of the kick velocity: 
130 <\V Z \< 220 km/s, and \V X \ < 115 km/sec. 

The parameter (3 is strictly less than one, reaching 
this limit only if no mass is lost in the second super- 
nova. With the conservative assumption that the presu- 
pernova mass of the companion star was at least 2M@, 

we have /3 ~ .8. If the misalignment of the system is 
S — 155°, then a very large kick in the direction oppo- 
site to the instantaneous orbital motion is required, with 
V y < —620 km/s for no mass loss and V y < —660 km/s 
for a 2Mq presupernova mass. The minimum total kick 
speeds are 630 km/s and 670 km/s. 

In the S ~ 25° case, much smaller kicks are acceptable. 
For the no-mass-loss case the constraint is — 80 < V y < 
11 km/s, with a minimum total kick speed of 149 km/s, 
for a 2Mq presupernova mass —120 < V y < —44 km/s 
with a minimum total of 173 km/s. 

Given the mass loss and natal kick, we can use equa- 
tions and |21 to calculate the systemic velocity kick ac- 
quired by the binary in the second supernova. In Fig. [21 
we show the relationship between the pre-supernova 
properties of the binary and the components of this kick 
perpendicular and parallel to the post-explosion orbital 
plane. It is immediately evident that a natal kick large 
enough to leave 8 ~ 155° inevitably leaves the binary 
system moving fast compared to the observed 122 km/s 
transverse velocity, suggesting that such large misalign- 
ments are only possible if the current binary velocity is 
fortuitously aligned very nearly along the line-of-sight 
to the pulsar. As we will see in the next section, such 
alignment is not consistent with our knowledge of the 
orientation of the binary. 

4.2. Motion of the binary in the Galaxy 
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Fig. 2. — The center-of-mass (or "systemic") velocity imparted 
to the PSR B1534+12 binary at the time of the second supernova, 
broken into components perpendicular and parallel to the post- 
explosion orbital plane. The sign of U± is unconstrained. The 
shaded regions mark the ini tial states allowed by the presupernova 
binary evolution models of Dcw i fc Polsl 1200.^ . The "island" of 
points at relatively low velocity are the solutions with S = 25°, 
while the island at high velocity are solutions with S = 155° . Rep- 
resentative error bars show the effect of varying <5 within the allowed 
la range. Lines of constant presupernova mass are shown, as are 
lines of constant presupernova orbital size, as parameterized by a 
(see eg. HI. 



The full three-dimensional orientation of the PSR 
B1534+12 orbit is known. The orbital inclination rel- 
ative to the plane of the sky (sini) is known from timing 
measurements. The usual ambiguity in the sign of cosi, 
which is not accessible from timing studies, is known 
from the profile observations. The angle of the line of 
nodes, corresponding to the rotation of the binary around 
the line of sight, is known from scintillation studies. To- 
gether, these observations fix the direction of the orbital 
angular momentum vector in space. Because tidal forces 
on the binary are very small, there has been no signif- 
icant change in the direction of the angular momentum 
vector since the current binary was formed in the second 
supernova explosion. 

The arguments of the last section give us, for any 
given presupernova parameters, the components of the 
systemic velocity kick parallel to and perpendicular to 
this angular momentum vector (e.g., Fig. [2J. Ideally, we 
could compare those kick velocity components directly 
to the observed space motion of the binary, but we must 
first overcome two significant problems. First, although 
the transverse motion of the binary is well determined, 
there is no direct constraint on the radial velocity. Sec- 
ond, the current velocity has been substantially affected 
by acceleration of the binary as it moves in the Galactic 
gravitational potential. 

Fortunately, the age of the binary is relatively low, so it 
is not unreasonable to follow its motion backwards to find 
possible birth locations and to thereby estimate the sys- 
temic kick at birth. It is very likely that the original pair 
of massive stars was born near the plane of the Galaxy, 
moving relatively slowly. The scale height of O type stars 
is small, ~ 50 pc (jStondll97^ . and their peculiar veloc- 
ities are low. Even after the first supernova explosion, 
the space velocity remains low, comparable to the frac- 
tional mass loss times the presupernova orbital velocity. 
In the canonical evolutionary model, where conservative 
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mass tran sfer prior to the first supernova leads to orbital 
widening ( Bha ttacharva fc van den Heuve]lll991I for ex- 
ample), a typical velocity kick is ~ 10 — 20 km/s. 

The velocity of the binary is therefore dominated by 
the recoil introduced in the second supernova explosion, 
when both the fractional mass loss and the pre-supernova 
orbital velocity are much higher than they were in the 
first explosion. The binary is currently far from the 
Galactic plane, at z = 680pc. Assuming the binary was 
moving with the local Galactic rotational velocity and 
was in the Galactic plane at the time of the second su- 
pernova will introduce only minor errors to estimates of 
the age and systemic kick velocity. 

We have studied the motion of th e pulsar through the 
Galax y using the potential model of iKuiiken &: Gilmord 
(1989). A series of trial radial velocities was used. In 
each case, the solar motion was added to the pulsar mo- 
tion to produce a space velocity. The orbit of the pulsar 
was then integrated backwards for 210 Myr. Each plane 
crossing represents a potential birth place and time for 
the binary system. Subtraction of the local rotational 
velocity then yields the systemic kick velocity needed at 
the time of the second supernova to carry the binary 
to its observed location. This can be decomposed into 
components parallel to and perpendicular to the post- 
supernova orbital plane, but we first discuss conclusions 
that are independent of the orientation of the system. 

Depending on its age, birth location, and kick velocity, 
we find that the binary could be on its first, second, 
third, fourth, or even (with careful fine tuning) even fifth 
excursion away from the Galactic plane. If this is not the 
initial departure from the plane, then we find that any 
pulsar age between 70 and 210 Myr is possible, as is any 
current radial velocity between about ±220 km/s. 

Radial velocities with v r ~ 200 km/s are only accept- 
able if it is still on its initial excursion away from the 
plane. A large velocity directed nearly radially away from 
us implies that the pulsar was born recently in the near 
neighborhood of the Sun. For example, v r = 200 km/s 
implies an age of just 7.9 Myr, while v r — 300 km/s 
implies an age of 4.5 Myr and a birth within a kilopar- 
sec from the Sun. While we cannot rule out arbitrarily 
large radial velocities, they require a fortuitous viewing 
angle and we believe them to be highly unlikely for other 
reasons. First, the young implied ages require that the 
pulsar birth period was very close to the observed pe- 
riod, and well below the spin-up line. This requires very 
fine tuning of the accretion parameters during spin-up. 
Second, unless we were extraordinarily lucky to find a 
double neutron star system born so recently in the solar 
neighborhood, the implied birthrate would be extraordi- 
narily high. 

If the binary has a large radial velocity towards us, 
then it must be returning to the plane for the first time. 
Again, fortuitous viewing alignments are difficult to rule 
out. However, in this case the velocity is directed gener- 
ally towards the Galactic center, so a radial velocity as 
high as —300 km/s implies a birthplace at a very large 
Galactocentric radius, R > 15 kpc. We regard birth at 
such large radii to be unlikely. 

We conclude, therefore, that the radial velocity is al- 
most certainly between —300 km/s and +220 km/s, and 
probably between ±220 km/s. We cannot set firm limits 
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Fig. 3.— Possible histories of the PSR B1534+12 binary. For 
each possible radial velocity, the orbit in the Galaxy was integrated 
backwards, and each Galactic plane crossing was recorded. The 
necessary systemic kick components parallel to and perpendicu- 
lar to the post-supernova orbit were calculated. Plotted solutions 
are those allowed for some choice of presupernova companion mass 
m,2i and semimajor axis a;, given 1<t ranges in the parameters 8 
and f2. Several families of acceptable solutions exist, marked by 
letters and summarized in Table HI Top left: Possible ages of the 
binary system, assuming the binary was born near the Galactic 
plane. The lowest curve in this figure (red points in the electronic 
edition) includes possible histories where the binary is leaving the 
plane for the first time (positive radial velocities) or returning for 
the first time (negative radial velocities). The second curve (green 
in the electronic edition) describes solutions where the binary is 
leaving or returning to the plane for the second time, and so forth. 
For the smallest allowed orbits, there is time in 210 Myr for the 
pulsar to be on its fifth excursion from the plane (having started its 
third complete orbit in the Galactic potential). Top right: Allowed 
presupernova companion masses. Bottom left: Allowed systemic 
kick velocities, relative to the postsupernova orbital plane. This 
figure can be compared with Fig. [5] Note that there are no solu- 
tions consistent with <5 155°. Bottom right: Magnitude of the 
asymmetric velocity kick imparted to the companion star at the 
time of the second supernova explosion. 



on the age of the system, but it is very likely more than 
about 100 Myr. Most important for our present pur- 
poses, we note that in all of the acceptable models the 
systemic kick given to the binary at birth was between 
~ 100 and ~ 240 km/s. 

Note that this conclusion is already strong enough to 
exclude the very large kicks needed to produce S ~ 155°. 
Also excluded are sce narios with large p resupernova 
masses, such as those of lDewi fc Polsl l|2003|) that avoid a 
second period of Roche lobe overflow, since the minimum 
presupernova mass in these models, ~ 4.5M©, requires a 
minimum systemic kick of ~ 270 km/s. The only models 
consistent with the observed low velocities are the case 
BB models in which the second helium star explodes with 
a mass of ~ 2 — 4M Q . 

4.3. Orientation of the orbit 

These constraints can be further tightened using our 
knowledge of the post-supernova orbital orientation. Af- 
ter propagating the binary back to a potential birth- 
place to find the direction and magnitude of the sys- 
temic kick that was needed at birth, we can resolve this 



Thorsett, Dewey, & Stairs 




-200 -1(10 '00 200 
V.^ (km/s) 



-200-100 100 200 
U M (km/a) 




100 ISO 200 
Up., (km/s) 



200 i no 100 200 
v ™a (km/s) 



Fig. 4. — Identical to Fig. |3 except 2<r ranges in S and £2 are 
accepted. The most important difference is that some second-pass 
solutions with larger negative radial velocity are allowed, somewhat 
broadening the range of allowed presupernova companion masses. 



TABLE 1 
Possible birth properties 



Key a Pass b 



Age c 
(Myr) 



m 2 i 



(km/s) 



A 


1 


40-60 


2.1- 


-2.7 


170- 


-260 


B 


1 


70-210 


3.0- 


-3.9 


230- 


-290 


C 


2 


70-210 


2.1 


-3.6 


170- 


-260 


D 


3 


130-190 


2.1 


-3.1 


170- 


-260 


E 


4 


170-210 


2.1 


-2.7 


170- 


-260 


F 


5 


200-210 


2.1 


-2.5 


170- 


-250 



a See Fig. 

b Excursion from the Galactic plane. For exam- 
ple, pass 2 means that the binary has already 
returned to the plane once since receiving its 
large velocity kick during the second supernova. 

c Time since second supernova explosion. 



kick into components perpendicular and parallel to the 
post-supernova orbital plane and solve equations and 
for the presupernova orbital size and companion mass. 
Using Kepler's law to eliminate a^, and some algebraic 
manipulation, we can rewrite equation [S] as a simple 
quadratic equation in rri2i'. 



= 



Gm\ s 



sin 2 8 



-2(m 2l + mi) 



(rri2f + m\)af 
ma/ + mi sin 2 8 



-m 2 f- 



sur 



8 Y G(m 2 f + mi) 



1 cos 8 



(14) 



Allowing for the sign ambiguity in u±, there are up to 
four real roots, each with a corresponding that can be 
determined from Kepler's equation. 

Physically interesting solutions have mn in the range 
2.1Mq-10Mq and ai within the range given by equa- 
tion 0] We consider 1 and 2a ranges in Q and 8 when 



computing the full set of possible progenitors. The re- 
sulting set of allowed solutions are shown in Figs.Eland^ 
We find viable solutions with 1-5 disk crossings and ages 
less than 210 Myr. There are no solutions with radial 
velocities less than —260 km/s or more than 240 km/s. 
The kick to the NS is relatively tightly constrained to 
between 170 and 290 km/s (between 150 km/s and 320 
allowing for 2a regions), while the progenitor mass can- 
not have been more than ~ 3.8 M (4.4M for 2a re- 
gions) , again implying that only case BB progenitors are 
possible. The kick direction is oriented within roughly 
20-40° or 140-160° of the pre-supernova orbital angu- 
lar momentum (with the ambiguity again arising from 
the ambiguity in the sign of the perpendicular compo- 
nent of the velocity kick) and the component of the natal 
kick that is in the presupernova orbital plane is directed 
roughly opposite to the instantaneous progenitor motion 
at the time of explosion. 

Several sources of uncertainty contribute to our results 
at a comparable level. Included in our analysis are the 
~ 20° uncertainty in f2, which introduces uncertainty in 
the decomposition of the birth velocity into parallel and 
perpendicular components, and the ~ 4° uncertainty in 8 
(e.g., Fig. [21). A smaller contribution is the 10 — 20 km/s 
uncertainty in the presupernova peculiar velocity of the 
binary. 

Another significant uncertainty, particularly for solu- 
tions with high space velocity, is the Galactic potential 
model used to model the binary's motion from the Galac- 
tic plane. We have repeated our analysis of the binary 
space motion using the three component potential of 
Uohnston eTall dl99^ . and find very similar results. The 
only significant difference is that for a very small range 
of radial velocities (between —205 and —220 km/s, on 
the second excursion from the Galactic plane, ages 200- 
210 Myr) there are models with more massive presuper- 
nova progenitors (77121 = 5-5.6Af Q ) that are consistent 
with the errors in f2 and 8. We regard these scenarios 
as unlikely, requiring careful fine-tuning as well as less- 
commonly used potential model, but the difference does 
highlight the systematic uncertainties in following the or- 
bit of such high-velocity objects in the Galaxy. (Those 
particular solutions are just returning to the plane from 
a maximum height of ~ 6 kpc.) 

In summary, a number of error sources as well as dis- 
crete ambiguities in the age of the system make it impos- 
sible to precisely determine the parameters of the presu- 
pernova orbit, even with the available data on the orbital 
orientation. However, the data are sufficient to strongly 
constrain the parameters: it is, for example, likely that 
the progenitor mass was in the range ~ 3 ± 1M , and 
that the natal kick is ~ 230 ± 60 km/s. 

5. CONCLUSIONS 

The combination of timing, precession, and proper mo- 
tion studies have allowed us to strongly constrain the evo- 
lution of the PSR B1534+12 binary. There are several 
conclusions that can be drawn: 

1. The misalignment of the pulsar spin axis from the 
orbital angular momentum is relatively small, with 
8 = 25.0 ± 3.8°. Although a misalignment of 
155.0±3.8° was formally allowed by observations of 
geodetic precession in the system, such a large mis- 
alignment would require an asymmetric kick at the 
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time of the second supernova explosion that is too 
large to reconcile with the current space velocity of 
the binary. 

2. Of the ev ol utiona ry scenarios studied by 
iDewi fc Polsl ((2003), only one is consistent 
with the observed properties and space velocity of 
the binary. The pulsar was the first-born neutron 
star of the pair. After an initial inspiral during 
either case B or C mass transfer as its companion 
evolved off the main sequence, the companion was 
a helium star between about 3.3 and 6.5M Q . A 
second phase of mass transfer (case BB) left the 
companion with mass 3 ± 1M© at the time of the 
second supernova explosion. Models that avoid a 
second phase of Roche-lobe-overflow-driven mass 
transfer in the binary by positing a high zero age 
helium main sequence mass for the companion 
also require an asymmetric kick that is too high. 

3. In the second supernova explosion, an asymmet- 
ric kick of 230 ± 60 km/s was required to produce 
the observed orbital elements. The direction of 
the kick relative to the orbital angular momentum 
vector (and presumably the angular momentum 
of the presupernova core) is tightly constrained, 
and is not obviously biased towards either the pole 
or the orbital plane. We note that the inferred 



asymmetric kick is somewhat larger than the 100- 
150 km/s most-likely kick f or PSR J0737-3039 
estimated from scintillation lRjmsoni i e^dJj200 " 
and e volutionary arguments ijWillems fe Kaloeer 
2004), and probably somewhat smaller than the 
kick given to the companion of PSR B1913+16, 
which is poorly constrained. An early analysis 
found most likely solut ions for PSR B19 13+16 to 
be ~ 300 - 500 km/s l|Wex et alJ 12000. prograde 
solution), but the range expands to 190-600 km/s 
when evolutionary sce narios including Roch e lobe 
overflow are included l|Willems et al.ll2004D . The 
limited measurement quality and complex selection 
biases make it impossible to estimate the distribu- 
tion of kick velocities imparted to newborn neutron 
stars in close binaries, but we note that there is no 
significant evidence in this small sample for a dif- 
ferent velocity distribution. 
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